Non-dynamic origin of the acoustic attenuation at high frequency in glasses 
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The sound attenuation in the THz region is studied down to T=16 K in glassy glycerol by inelastic 
x-ray scattering. At striking variance with the decrease found below ~100 K in the GHz data, the 
attenuation in the THz range does not show any T dependence. This result i) indicates the presence 
of two different attenuation mechanisms, active respectively in the high and low frequency limits; ii) 
demonstrates the non-dynamic origin of the attenuation of THz sound waves, and confirms a similar 
conclusion obtained in Si02 glass by molecular dynamics; and Hi) supports the low frequency 
attenuation mechanism proposed by Fabian and Allen (Phys. Rev. Lett. 82, 1478 (1999) ). 



PACS numbers : 78.35.-l-c, 78.70.Ck, 61.20.Lc, 83.50.Fc 



One of the most important and still unsettled subjects 
in the physics of topologically disordered systems regards 
the mechanisms for the propagation and attenuation of 
density fluctuations. The propagating nature of acoustic 
waves, as seen by Ultrasonic and Brillouin Light Scatter- 
ing (BLS) measurements in the MHz and GHz region 
respectively, has been shov^rn to persist up to the THz 
region by the existence of a linear relation between the 
peak energy, E, and the momentum transfer, Q, of the 
inelastic features observed in the dynamic structure fac- 
tor, S{Q, E), of glasses 0]. This result is the outcome of 
extensive studies on the shape of S{Q, E) performed us- 
ing Molecular Dynamics (MD) simulations j^-^ and the 
newly developed Inelastic X-rays Scattering (IXS) tech- 
nique This latter technique allows to study the 
S(Q,E) in the "high" Q range (Q «1-10 nm"!), thus 
increasing by about two orders of magnitude the Q val- 
ues typically investigated by BLS [Q «0. 01-0. 04 nm~^). 
In the IXS and MD Q-range, beside the persistence of 
a linear dispersion of the acoustic excitation energies, 
one also observes a progressive broadening of the inelas- 
tic features, which is responsible for their disappearance 
at a certain Qm value. Typically Qm is some tenths of 
Qm — the position of the first sharp diffraction peak in 
the static structure factor, S{Q) [Q. The study of the 
mechanisms leading to this damping, and, therefore, the 
investigation of the sound waves attenuation at these Q- 
values — characteristic of structural correlations at the 
interparticlc level — is obviously of great interest. 

The acoustic excitations at frequencies in the THz 
range, as measured so far in glasses and glass forming liq- 
uids by IXS, have a linewidth parameter Fq which seems 
to show a dependence [Q . Moreover — in all the IXS 
data reported so far — Tq/Q^ has a negligible temper- 
ature dependence in a wide temperature region ranging 
from values well below the glass transition temperature, 
Tg, up to the liquid phase [|l|. At variance with this be- 
havior, as well known, the linewidth of the excitations in 



the GHz region, measured by BLS, show a relevant tem- 
perature dependence, which becomes particularly strong 
in the limit of very small temperatures |p|-p^. The tem- 
perature dependence of the linewidth in the GHz range 
has motivated many theoretical studies, leading to differ- 
ent hypotheses on the frequency (or Q) evolution of the 
attenuation mechanisms p^ ]. 

In this Letter we report an IXS study on the low tem- 
perature behavior of the excitations linewidth in glassy 
glycerol. Specifically, we concentrate on the study of 
THz excitations in the temperature region where the 
BLS data in the GHz range show a marked tempera- 
ture variation. Within the error bar, the linewidth mea- 
sured by IXS is temperature-independent in the whole 
O.lTg to Tg region, whereas, in this same region, the BLS 
linewidth increases by more than a factor of ten. This 
two opposit behaviors indicate that there are at least two 
different attenuation mechanisms: i) One of dynamic ori- 
gin dominant in the low Q (low frequency) region, and 
it) A second one, dominant at high Q, whose tempera- 
ture independence suggests that its origin is due to the 
structural disorder of the glass. The glycerol results are 
confirmed by a similar sound attenuation behavior found 
in vitreous silica, as obtained by the analysis of existing 
BLS, IXS and MD data. The observation of two distinct 
attenuation mechanisms, each one dominant in a different 
Q region, implies the existence of a cross-over frequency, 
which lies in the 100 GHz range for both of the studied 
glasses. It also suggests that the frequency dependence of 
the dynamic contribution to the sound attenuation agrees 
with the one recently predicted by Fabian and Allen . 

The experiment has been carried out at the new very 
high energy resolution IXS beamline ID28, at the Euro- 
pean Synchrotron Radiation Facility. The incident x-ray 
beam is obtained by a back-scattering monochromator 
operating at the Si(ll 11 11) reflection ||l^. The scat- 
tered photons are collected by a spherical silicon crystal 
analyzer, also operating at the Si(ll 11 11) reflection [fsl. 
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The monochromatic beam has an energy of w 21,748 eV 
and an intensity of 2-10® photons/s. The total energy res- 
olution — obtained from the measurement of S{Qm,E) 
in a Plexiglas sample which is dominated by elastic scat- 
tering — is 1.5 meV full- width-half-maximum (fwhm). 
The momentum transfer, Q — 2kosin{9s/2), with ko the 
wavevector of the incident photon and 0^ the scattering 
angle, is selected between 2 and 4 nm~^ by rotating a 7 m 
long analyser arm in the horizontal scattering plane. The 
total Q resolution has been set to 0.2 nm~^ . Energy scans 
are done by varying the relative temperature between the 
monochromator and analyzer crystals. Each scan took 
about 180 min, and each (Q, T)-point spectrum has been 
obtained from the average of 2 to 8 scans depending on 
the sample temperature. The data have been normalized 
to the intensity of the incident beam. The sample cell is 
made out of a pyrex-glass tube (4 (10) mm inner (outer) 
diameter and 20 mm length), capped with two diamond 
single crystals discs, 1 mm thick, to minimize undesired 
scattering signals. The cell has been loaded with high pu- 
rity glycerol in an argon glove box. In the Q — E region of 
interest, empty cell measurements gave the flat electronic 
detector background of 0.6 counts/min. The cell length 
was chosen to be comparable to the x-ray photoabsorp- 
tion length, and multiple scattering was negligible. 

The spectra have been collected at T=16, 45, 75, 114, 
145 and 167 K, and, as examples, those at T = 16 
and 167 K are reported in Fig. 1 for different Q-values. 
The full lines are the fits to the data, obtained us- 
ing a model function made by the convolution of the 
experimentally determined resolution function with a 
delta function for the elastic peak and a Damped Har- 
monic Oscillator (DHO) model for the inelastic peaks 
[01 . This model for the S{Q,E) results from the as- 
sumption that the memory function, mQ{t) fl^ , en- 
tering in the Langevin equation for the considered Q- 
component of the density fluctuation, has a time de- 
pendence as: mQ{t) = 2TQ{T)6{t) + A^(T) The 
presence of the structural a-relaxation, observed in the 
liquid state and frozen in the glass, and of other relax- 
ation processes with characteristic times slower than « 1 
ps, is reflected in the parameter Aq(T), whose value de- 
termines the change of the sound velocity, c, between 
the fully relaxed (co) and unrelaxed (coo) limiting val- 
ues: A^(T) ^ Q^{cl^ - cl). The parameter Tq{T) de- 
termines the width of the side peaks, i. e. the sound 
wave attenuation coefficient, a = 2'KTQ/hc. However, 
fits made with the DHO model or with different fitting 
function gave values for the FWHM of the inelastic peaks 
consistent among each other within their statistical un- 
certainties, indicating, therefore, the insensibility of the 
results to the specific model for the inelastic peaks. As 
it is evident already from the raw data, Tq{T) has a 
marked Q-dependence while its T-dependence, if any, is 
much smaller. This is better seen by the dotted lines in 
Fig. 1, which represent the unconvoluted inelastic part of 



S{Q,E). 

The values of Fg (T) resulting from the fit of the IXS 
data of Fig. 1 are reported as a function of Q in Fig. 2. 
In the same figure are also shown the Tq{T) obtained 
from IXS measurements at 175 K [Q, and those obtained 
from literature BLS spectra measured at Q «0.03 nm~^ 
and temperatures similar to the IXS' ones Ppl|]. This 
figure demonstrates that, within the error bars, Tq{T) is 
T-independent in the Q region covered by IXS. On the 
contrary, Tq (T) shows a marked T-dependence at the Q 
value of the BLS measurements. These T-dependencies 
are emphasized in Fig. 3, where Tq{T)/Q'^ is plotted as a 
fimction of T for Q in the BLS region (Q « 0.035 nm~^) 
§,|ll and for Q in the IXS region. Here the crossed sym- 
bols refer to IXS measurements at fixed Q {Q = 2 nm~^) 
and the full symbols to the average of Vq{T)/Q^ over 
the (5 = 2 — 4 nm~^ region. This figure confirms that, at 
the high Q values, Vq{T) is substantially constant even 
at very low T, whereas, at the low Q values, it increases 
with temperature up to k, 100 K, where it seems to reach 
a plateau. The further increase above Tg is due to the 
a-relaxation; it is only seen in the low Q data as it would 
affect the high Q data at higher T. 

The specific dynamic mechanisms (anharmonicity, re- 
laxation processes, floppy modes, two level systems...) 
at the origin of the acoustic attenuation observed at the 
BLS' Q values, as well as their temperature dependence, 
have been widely investigated in the past jl5|. In con- 
trast to what it is found in the BLS' Q region, the be- 
havior of Tq{T) in the high Q region, as reported in 
Figs. 2 and 3, shows that here the sound attenuation is 
not determined by temperature activated dynamic pro- 
cesses. Consequently, in this Q range, at variance with 
the crystalline state where the absence of dynamic pro- 
cesses would imply no sound attenuation, in the glass the 
observed non vanishing value of Fg must have a " struc- 
tural" origin, i. e. it must be due to the topological dis- 
order of the glass structure. 

The picture coming from the reported data suggests 
that one can express rg(r) = F*^' (T) + F<|', where 
Fg'' (T) is a temperature dependent dynamic part and 
Fg' is due to topological disorder. The Q-dependence of 
Fg'' (T) must be such to be the dominant term of Fg(r) 
at small Q, while it must be negligible at large Q. This 
behavior is consistent with a recent calculation of the dy- 
namic (anharmonicity) contribution to the sound atten- 
uation by Fabian and Allen , who predict a Fg (T) oc 
up to a Qc value, above which Fg(T) = const. In the 
case of amorphous silicon, Qc has been calculated to be 
in the 0.1 nm~^ range [|l^. Although there is no such 
calculation for glycerol, the present results indicate that 
also in this glass the crossover takes place at Q values 
between the BLS' and IXS' ones, i. e. in the Vc «100 
GHz frequency range. It is worth to note that a relax- 
ation process with characteristic time r = 1/2'kVc ~ 2 ps, 
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and responsible for a linewidth of the order Tq/Q"^ « 0.2 
meV/nm~^, should also give a dispersion of the sound 
velocity, 5c, given by 5c/ c — TrVq/Q'^Thc^ ~ 1%, a value 
too small to be detectable with the accuracy achievable 
at present. Therefore, in the present case, the change of 
sound velocity cannot be used to estimate the crossover 
value Qc- 

The Q-dependence of Fg ' , as already observed before 
in many other glasses and glass forming systems [Q, is 
well represented by a law, Fq' = D Q^, in the Q 
region covered by IXS. This law, shown as a full line 
in Fig. 2, however, cannot be extrapolated (thin full line) 
to low Q values because it would predict width values in 
excess to the measured Fg (T) . This observation excludes 
that F^' oc Q2 in the whole 0.01-10 nm"! Q range. Un- 
der the hypothesis that Fg' = Q^, and assuming that at 
the lowest measured temperature Fg' = Fg, one finds 
consistency within the error bars of both the IXS data 
and the BLS low temperature point, as shown by the 
dashed line in Fig. 2 obtained with 7 = 2.5. This esti- 
mate of 7 is a low limiting value, because the BLS low T 
width could still be partially affected by a dynamic con- 
tribution. It is not clear, however, whether the hypoth- 
esized power law indeed provides a good representation 
of Fg ' , some hints on this issue can be gathered by the 
study of other glasses. 

The previous picture for the sound attenuation in the 
glycerol glass is further substantiated by the existing IXS 
|2|-|§, BLS ]T0|-P, p5|, M P and Picosecond Opti- 
cal Technique (POT)pq], data in another prototypical 
glass: vitreous silica {v-Si02)- We report in Figs. 4 and 
5 the Fg(T) values for Y-Si02 in a format equivalent to 
that of Figs. 2 and 3. In Fig. 4, the Fg(T) values in the 
1 to 5 nm~^ region, as obtained by IXS, do not show 
any relevant T-dependence in the 300-1450 K range. As 
a consequence of contrast problems due to the limited 
energy resolution of the IXS spectrometer, at tempera- 
tures below 300 K, it has not been possible to discrimi- 
nate the inelastic signal from the tails of the elastic one. 
The missing low temperature IXS data are supplied by 
an extended MD simulation performed in the harmonic 
approximation (T = K), which provides a S{Q, E) line- 
shape in excellent agreement with the IXS data, and con- 
firms the absence of any relevant T-dependence of the 
sound attenuation in the whole T=0-1450 K and (5=1- 
5 nm^^ ranges. On the contrary, as in glycerol and as 
emphasized in Fig. 5, the BLS data of v-Si02 data show 
a large T-dependence. As in other glasses, also the high 
Q v-Si02 data of Fg have a Q-dependence well repre- 
sented by a law (full line in Fig. 4). In \-SiO2 are 
also available room T data in a Q region intermediate to 
BLS and IXS (Q=0.03-0.4 nm^i), as obtained by POT. 
These data, however, show a relevant inconsistency with 
BLS data measured in the same Q and T ranges. 

The similar behavior between Y-Si02 and glycerol al- 



lows to formulate also for \-SiO2 the same hypotheses 
on the Q and T-dependencies of Fj^' (T) and F<^' . It is 
worth to note, however, that the ability of the power law 
for Fg' (which should have 7 — 2.6) to pass through all 
the IXS data and the low T BLS point is substantially 
worse than in glycerol. Therefore the y-Si02 results 
clearly indicates that the power law hypothesis is wrong 
and that Fg' has a more complex behavior, namely it is 
oc in the high Q regime (Q > 1 nm~^) and it has a 
steeper behavior at low Q values. 

In conclusion, we have shown that in glycerol and sil- 
ica glasses the dominant sound attenuation mechanism 
has a different origin in the Q region spanned by the 
IXS or BLS techniques. The temperature independence 
of the attenuation at large Q, corresponding to frequen- 
cies in the THz region, implies that its origin is struc- 
tural and it is due to the disorder. On the contrary, the 
well known strong T-dependence found in the low Q re- 
gion, at GHz frequencies and below, implies a dynamic 
origin of the sound attenuation ||l^] . These findings im- 
ply a cross-over region between the two regimes, which 
should lie in the 100 GHz frequency range. The overall 
T and Q-dependencies of the attenuation considered here 
is consistent with a dynamic (F^* ) part that closely fol- 
lows the one proposed by Fabian and Allen jl^, and a 
structural part (Fg ' ) that has a behavior at high Q 
{Q > 1 nm^^) and a steeper Q dependence at lower Q. 
The hypothesis on a cross-over between the two different 
attenuation mechanisms considered here call for further 
studies, where the Q and T-dependencies of the sound 
waves is thoroughly investigated using the BLS and IXS 
technique in a wider number of glass materials. 
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FIGURE CAPTIONS 

FIG. 1 - Inelastic x-ray scattering spectrum of glycerol at T=16 
and 167 K and at the indicated Q-values. The full lines are 
the best fit to the data as discussed in the text. The dotted 
lines represent the uncovoluted inelastic contributions to the 
fit. 

FIG. 2 - The lincwidth parameters Fq for glycerol are reported 
as a function of Q at the indicated temperatures in the IXS' 
(full symbols) and BLS' (open symbols) Q regions. The 
inset shows an enlargement of the IXS' Q region. The full 
line represents the behavior, which is the best fit to Tq 
at high Q. Also shown in the figure axe extrapolation of 
the law in the low Q range (thin full line) and the Q^'^ 
dependence (dashed line) indicated by the low T BLS data. 
The inset shows an enlargement of the high Q region. 

FIG. 3 - Temperature dependence of Tq/Q^ in glycerol at Q as 
0.03 nm^-'^ (open symbols), at Q = 2 nm^-'^ (crossed sym- 
bols) and averaged over the (5 = 2 — 4 mn^^ region (full 
symbols). The vertical dashed line indicates the glass tran- 
sition temperature, rg=187 K. 

FIG. 4 - Same as in Fig. 2 but for vitreous silica. Data from MD 
simulation (crossed symbols) and from POT (stars) arc also 
reported. Here the low T BLS data indicate (if any) a power 
law with exponent 7 = 2.6. The inset shows an enlargement 
of the high Q region. 

FIG. 5 - Same as in Fig. 2 but for vitreous silica. The inset 
shows an enlargement of the low temperature region. The 
open circle refer to BLS data at Q Ri 0.035 (□ and o) and 
Q ss 0.025 (O) nm-i, the full squares to IXS data at Q = 
1.6 nm~^, and the crossed circle to MD simulation in the 
harmonic (T= K) limit. The inset shows an enlargement 
of the low T region. 
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